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Abstract. We give a one-to-one correspondence between classes 
of density matrices under local unitary invariance and the dou- 
ble cosets of unitary groups. We show that the interrelationship 
among classes of local unitary equivalent multi-partite mixed states 
is independent from the actual values of the eigenvalues and only 
depends on the multiplicities of the eigenvalues. The interpretation 
in terms of homogeneous spaces of unitary groups is also discussed. 
PACS numbers: 03.67.Mn, 02.20.Hj, 02.10.Ox 



1. Introduction 

Theoretic foundation of quantum computation and quantum infor- 
mation [NC] has been postulated through the language of quantum 
mechanics. At the beginning of the development of quantum mechan- 
ics Schroginger pointed out the unitary freedom of the quantum system 
[SI in], which is one of the spotlights in the whole theory and has been 
used in many applications in later developments. In |HJWj the free- 
dom of purification of quantum states was studied with help of Schmidt 
decomposition. Furthermore the probability distributions given by a 
given density matrix are characterized by majorization [N]. These re- 
sults have played an important role in quantum statistics, quantum 
computation and quantum information. 

Quantum entanglement is one of the key issues in quantum computa- 
tion, and local unitary invariance is an importance aspect in quantum 
entanglement and their applications. For instance, several fast algo- 
rithms are discovered based on special properties of local invariance 
and quantum entanglements. The maximum entanglement can also 
be achieved via local unitary action. Local unitary transformations 
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have also been used for fractional entanglement and in other prop- 
erties (cf. |HHHj ). In various investigations of local unitary equiv- 
alence lots of efforts have been made in seeking possible invariants 
[XFGl [run [ZSBl lACFWj . There have also been results on non-local 
unitary equivalence [DCj in exploring feasible quantum gates. In |YLF] 
a geometric class was proposed to study the problem for two partite 
quantum states and it has been shown that two quantum states are 
equivalent if and only if their representation classes are the same. In 
[FJj we have proposed a new operational method to study local unitary 
equivalence of density matrices and the relationship with separability. 

It is clear that in classification of equivalent density matrices one first 
needs to see whether two density operators have the same eigenvalues 
or not. Therefore one should separate the easy task of comparing 
eigenvalues from the problem and focus on other invariants of unitary 
operations. In the current paper we apply this strategy to a larger 
quantum system consisting of several quantum sub-systems and study 
their unitary properties from the viewpoints of both global and local 
pictures. 

We first look at density matrices with specified spectrum of eigen- 
values and study the properties of the quantum systems using group- 
theoretic methods. The local unitary properties are investigated with 
the help of fixed point subgroups of the concerned Lie group and we 
then set up a one-to-one correspondence between classes of local unitary 
equivalences and double cosets of the unitary group U{n) by the ad- 
ditive and multiplicative Young subgroups. Next we recast the space 
of the equivalence classes in terms of induced representations of the 
Young subgroups and show their deep connections with homogeneous 
spaces of unitary groups. 

One of our main results shows that local unitary equivalence depends 
only on the multiplicities of the eigenvalues and is independent from the 
actual values of the eigenvalues. This means that with given spectrum 
of the density matrices the local unitary equivalence is determined by 
the combinatorics of the quantum states. We hope that the combina- 
torics and invariant theory will shed light in further investigation of 
local unitary properties and entanglement of quantum states. 



2. Local unitary equivalence 

Let if be a Hilbert space affording the quantum system. If the 
quantum system is in a number of normalized quantum states 
with probabihty Pi, then the density operator for the quantum system 
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is the convex combination 

P = ^Pi\'^i){'ilJi\, 

i 

where Y2iPi — 1- Equivalently density operators are characterized as 
non-negative operators on the Hilbert space with unit trace. 

Among all physically possible quantum systems afforded by H, one 
needs to judge whether two systems are distinct or not. To study 
the equivalence of various quantum systems on H, one only needs to 
consider the unitary linear group U{n) over C. Suppose two density 
matrices pi have the same set of eigenvalues and their sets of multi- 
plicities for each eigenvalue arc also the same. Let the corresponding 
(orthonormal) eigenstates be \(j){X'f^)) and |-?/'(Aj-''^)) respectively, where 
j — 1, ■ ■ ■ ,mi, the multiplicity of the eigenvalue Aj. Then we can write 

Pi = Y.\m^^)){<i^i^\ 

P2 = $:a#(a?'))(V^(a?))|. 

ij 

We define g to be the hnear transformation on H sending \(f)(X^^)) i-^ 

\ijj{Xl^^^)), then g G U{n) due to orthogonahty of the eigenstates. Then 
one easily checks that 

P2 = gpig"^ = gpig^- 

In other words density matrices with similar spectral decomposition 
are unitary equivalent and all physical features of the quantum system 
are captured by properties of the unitary group U (n) . 

The above discussion can also be viewed in terms of group action. Let 
D{n) be the set of density matrices on the Hilbert space H of dimension 
n. The evolution of the density operator p is given by p — > gpg\ 
where g = g{t) e U{n) and f is transpose and conjugation. Clearly 
gpg^ is positive and tr{gpg^) = tr{p) = 1. We will adopt the standard 
convention in group theory to write hgh~^ — g^ for g,h & G. 

Lemma 2.1. The kernel of the action U{n) x D[n) — > D[n) is 
{e'^I\e e R} < U{n). 

Proof. Suppose gpg'^ = p for all density matrices p. If p is any diagonal 
matrix with distinct eigenvalues, then g must be diagonal by compu- 
tation. A generalized (or signed) permutation matrix is one that each 
row or column has only one and only non-zero entry and this non-zero 
entry is 1 or —1. One notes that any matrix commutes with a gen- 
eralized permutation matrix must be diagonal. To show that g — al 
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with a = e a complex number of unit modulus, we can take positive 
generalized permutation matrix, then all diagonal entries of g must be 
equal. □ 

For p G D{n) we define the invariant subgroup Up by Up = {g E 
U{n) \ gpg'^ = p}. As g E U{n) is unitary: gg"^ = g'^g = I, we have im- 
mediately that gpg^ = p iS g p = pg . The next statement is immediate 
from definition. 

Lemma 2.2. For any h E GL{n) one has Uph = hUph^ = iUp)^ . 

Let p be a density matrix of size n x n and thus is diagonalizable 
due to hermiticity. According to basic theory in linear algebra [L] there 
exists a unitary matrix g E U{n) such that 

p = gdiagiaJx-,, ■ ■ ■ ,ailxi)g\ 

where Aj are non-negative integers such that Ai + ■ ■ ■ + A; = n. We can 
further suppose that A, are arranged in a descending order: Ai > ■ ■ ■ > 
Xi > 0, thus A is a partition of n. Here / is the number of parts of A. 

Definition 2.3. We say a density matrix p is of type A, a partition of 
n, if p has eigenvalue multiplicities: Ai, ■ ■ ■ , A/, where Ai H — ■ + A/ = n. 

Two density matrices are equivalent if they are related by a unitary 
transformation: pi = gp2g\ g E U{n). This is in agreement with the 
equivalence relation given by the group action of U{n) on D{n). Thus 
the invariant subgroups of equivalent density matrices are conjugate in 
U (n) . We will call the set of equivalent density matrices of type A the 
equivalent class of type A, and denote the class by [A]. 

Proposition 2.1. If p = diag{ailx^, ■ ■ ■ i^ilxi) is a density matrix, 
then Up = U{Xi) x • • ■ x U{Xi). More generally if p is of type X = 
(Ai, ■ ■ ■ , A;), then Up = gU{Xi) x ■ ■ ■ x U{Xi)g^ for some g E U{n). 

Proof. If p = diag{ail\^^ ■ ■ ■ , ailxj^), it follows from direct computation 
that Up = U{Xi) X ■ • • X U{Xi). Then for a general density matrix p 
of type A, we have that p = gdiag{ailx^, ■ ■ ■ , ailxi)g'^ {g E U = U{n)), 
then the result follows from Lemma [2. 2 [ □ 

The classes of density matrices can be characterized by the group 
action. 

Theorem 2.4. Let [A] be the equivalent class of density matrices of 
type X, then [A] ~ U/Up, where p is some density matrix of type X. 
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Proof. Let U = U{n). The action U x [A] — > [A] given by {g, p) ^ gpg'^ 
is transitive. As the class [A] = {p^\g ^ U} for a fixed density matrix 
p, we see that 

therefore U/Up ~ [A(p)]. It follows from Proposition 12.11 that U/Up = 
U/gU^p)g-^ ^ □ 

We now consider local unitary equivalence. Let Hi be two Hilbert 
spaces with dimensions rii. We say two density operators pi on the 
space H = Hi® H2 are equivalent under local transformation iff pi = 
{91 ® g2)P2{.gi ® 92)"^ for some gi G Uirii) < End{Hi). We recall that 
are (globally) equivalent if pi = gp2g'' for g G U{n) < End{H), where 
n = nin2. More generally, for multi-partite cases the local unitary 
group is Un = U{ni) ® ■ ■ ■ ® U{nr), where Ui = dim{Hi). 

Theorem 2.5. Let H = Hi ® ■ ■ ■ ® H^ be the global Hilbert space 
with dim{Hi) = rii, and as above. Then the type A equivalent multi- 
partite mixed states under local equivalence are in one-to-one correspon- 
dent to double cosets of the unitary group U{n) by the Young subgroups: 
Un\U/Ux, where Ux = U{^i) x ■ ■ ■ x f/(A;). Moreover, the representa- 
tive of the double coset determined by the density matrix p is given by 
the matrix of the orthonormal eigenstates of p. 

Proof. Globally equivalent density matrices are determined by their 
eigenvalue spectrum and multiplicities. The class of (locally) equivalent 
density matrices of type A consists of 

where A = diag{ail\^, ■ ■ ■ , ail\^) and g \s a. unitary matrix consisting 
of orthonormal eigenstates of p. Here A is fixed. When two density 
matrices pi = giAgj are local unitary equivalent, then pi = kp2k~^ for 
some k & Un- It then follows that 

giAgl = kg2Aglk^ =^ 9i^kg2 = c E Ua, 

and Ua = U\ = U{\i) x ■■■ x U{\i) by Proposition 12.11 Therefore 
UngiUx = Ur,g2U\. 

Conversely, suppose two double cosets UngiU\ are the same. Then 
gi = kg2C for some k eU^ and c E Ux. It follows that gi^kg2 = G 
Ux, thus giAgl = kg2Ag\k^ . It is clear that pi = giAg\ gives rise to 
a density matrix and p2 = g2Ag\ gives rise to a globally equivalent 
density matrix (as pi = p^. Their corresponding classes (under local 
equivalence) are also equal: 

[^7iA^I] = [(72A^^]. 
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□ 



3. Induced representations 



Let A be a partition of n, n = Ai + ■ ■ ■ + A^, and n be a factorization 
oi n, n = rii ■ ■ - fir, where n is the dimension of the underlying (global) 
Hilbert space. The Young subgroup Ux associated with the partition A 
is the direct product f/(Ai) x ■ ■ ■ x U{\i). From our previous discussion 
it is clear that one also needs to consider the multiplicative Young 
subgroup Un = U{ni) ® ■ ■ ■ ® Uirir) associated to the factorization n 
of n. The additive Young subgroup U\ = U{\i) x ■ ■ ■ x f/(A;) can be 
imbedded into U{n) in the canonical manner: 



while the multiplicative Young subgroup Un = U{ni) ® ■ ■ ■ ® U{nr) is 
imbedded into U (n) via tensor product 



As we have seen in Section [2] the double cosets given by the Young 
subgroups U\ and multiplicative Young subgroup f/n are in one-to-one 
correspondence of quantum multi-partite systems with dimensions n^. 
We now reformulate this correspondence in terms of representations of 
the unitary group. We can restrict ourself to the special unitary group 
without loss of generality. Let G = SU{n) and H = SU\, and it is 
clear that H is a. closed subgroup of the compact Lie group G. 

Let C be the trivial representation of the subgroup H. The induced 
representation IndfjC, as a vector space, is the space of all complex 
continuous functions on H satisfying the condition: 



The induced representation Ind'^C is isomorphic to C°(G'/if, C), the 
space of continuous functions on the cosets G/H |BDj . Here the action 
is the left translation: L{g)f{x) = f{g^^x). The canonical basis of 
the representation is given by characteristic functions (pgH, indexed by 
right cosets G/H. Here 




gi ® ■ ■ ■ <S) gr ^ Un ^ gi ^ ■ ■ ■ ® gr & U{n). 



f{gh) = h-'f{g), geG,heH. 
The action of G on the induced representation is given by 
(3.1) ig-f)ix) = f{g-'x), g,xeG. 



(3.2) 
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The induced representation Ind'§C is essentially equivalent to the set 
of unitary classes of density matrices as seen from the following obser- 
vation. 

Example 3.1. Suppose n = 1 ■ n is the factorization of n, then local 
unitary equivalence coincides with global unitary equivalence. The dou- 
ble cosets U-n\U/U\ reduces to cosets U/U\. The unitary equivalence 
classes (with same set of spectrum) are also in one-to-one correspon- 
dence of the set of partitions Vin): 

n = \l^ h A;. 

Each equivalence class is represented by a homogeneous space of the 
unitary group U{n) determined by the partition A. The generating 
function of the number of unitary equivalence classes is given by [Mj 

oo °° 1 

^ #(equiv class) g" = JJ ^_ ^ - 

n=0 n=l ^ 

Example 3.2. Density matrices on any Hilbert space of prime dimen- 
sion are always inseparable. Multipartite states should live on Hilbert 
spaces of dimension equal to powers of prime numbers. 

Local unitary equivalent classes can also be studied from represen- 
tation theoretic viewpoints. Let K = SUn corresponding to the factor- 
ization n, and H = SUx < G = SU{n) as above. Both H and K are 
closed subgroups of G. We consider the restriction ResxInd'^C Let 
W be any if-module, then it follows from Frobenius reciprocity |BD] 
that 

HomK{ResKlnd%C, W) ~ Homcilnd^C, Ind%W) 

~ HorriHiC, Resnlnd^W). 

In particular, when W = C, one has the duality ResxInd'^C ^ 
Resnlnd'^C The linear operators in HomclInd'fjC, Ind^C) in gen- 
eral are certain distributions on the direct product of the unitary group 
by Schwartz's distribution theory [Scj. If the intertwining number (the 
dimension of Homcilnd^C, Ind^W)) is finite, then they are equal to 
the number of double cosets of K\G/H by Mackey's theorem [Mcj . 
Thus we obtain the following result. 

Theorem 3.3. Let H = Hi® ■ ■ ■ ® he the global Hilbert space with 
dim{Hi) = Hi, and as above. Then the classes of globally equivalent 
multi-partite mixed states under local equivalence are isomorphic to the 
restriction of the induced representation Ind^^'^^^'Q to the subgroup SUn- 
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4. Homogeneous spaces 

Homogeneous spaces are important non-empty manifolds with a tran- 
sitive action of a Lie group. The quotient space G/H are special ex- 
amples of homogeneous spaces. If one views equivalent local unitary 
classes as points in the homogeneous spaces one may understand the 
geometric meaning of the statement that all points are the same. 

Let n be a fixed prime decomposition of n and A be a fixed partition 
of n. If two density matrices pi are locally equivalent, then their cor- 
responding double cosets U^giUx and U^giUx are the same. Note that 
in this description the equivalence relations are determined completely 
by the double cosets, and no information about the actual eigenvalues 
are needed. In other words, classes of density matrices with the same 
eigenvalue distribution can be mapped to classes of density matrices 
with the same eigenvalue multiplicities. Therefore we have proved the 
following result. 

Theorem 4.1. Let H = Hi® ■ ■ ■ ® H^. be the global Hilbert space with 
dim{Hi) = Hi, and Un as above. Then the class of globally equiva- 
lent multi-partite mixed states under local equivalence can be mapped 
to another class of globally equivalent multi-partite mixed states with 
the same eigenvalue distribution and multiplicities. The two classes 
are both described by the homogeneous spaces Un\U/Ux, where the par- 
tition A is given by the eigenvalue multiplicities. 

This result shows that the local unitary equivalence relations does 
not depend on the actual values of the eigenvalues, and all classes can 
be viewed as "equal" in geometric sense. Therefore we can parameterize 
the classes of density matrices by the set 

X {U^\Uin)/Ux), 

where M is used for the eigenvalue. 

Example 4.2. For e > and 0</<l — ewe consider the two qubit 
Werner state 



P 



1+2/ 1-4/ 



6 6 



When e = 0, this is the usual Werner state |Wj. The eigenvalues are 

(l-/ + e)/3,(l-/)/3,(l-/-e)/3,/. 

The nontrivial factorization is certainly 4 = 2-2 and the partition of 
the eigenvalue multiplicities are A = (1111) for e > and / ^ ^j^. 
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When e > and / = or 1/4 then A = (211). If e = 0, then 
A = (31) for / 7^ 1/4 and A = (4) for / = 1/4. From the picture there 
are essentially one dense class of local unitary equivalence; two one- 
dimensional classes of unitary equivalence; and one degenerate class 
of unitary equivalence. The most interesting cases are A = (31) and 
(21^), where the local unitary classes are classified by double cosets of 
SU^22)\SUiA)/{SUi3) X SU{1)) and 5f/(22)\5f/(4)/(5f/(2) x 5f/(l) x 
SU{1)) respectively. All classes with the same partition are viewed as 
equal. 



/ 




5. Conclusion 

We have established a one-to-one correspondence between the set of 
density matrices of multi-partite systems and the set of double cosets of 
additive and multiplicative Young subgroups of unitary groups. Our re- 
sults show that the interrelationship among the classes of local unitary 
equivalent multi-partite mixed states is independent from the actual 
values of the eigenvalues and only depends on the multiplicities of the 
eigenvalues. This interesting phenomenon allows us to look at the local 
unitary equivalence for all global unitary density matrices at the same 
time. The geometry of the classes are given by that of homogeneous 
spaces of the unitary group as well as the invariant theory of classi- 
cal groups (cf. [HJ). It is expected that the combinatorics of density 
matrices and invariant theory will also play a role in this study. 
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